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Abstract. The Lunar Atmosphere Dust Environment Explorer (LADEE) spacecraft is being designed for a mission 
featuring low altitude orbits of the Moon to take relevant ambient measurements before that environment becomes altered 
by future exploration activities. Instruments include a neutral mass spectrometer capable of measuring ambient species 
density levels below 100 molecules/cm 3 . Coincidentally, with a favorable combination of spacecraft orientations, it is 
also possible to measure plume gases from LADEE attitude control system thruster operations as they are reflected from 
the daytime lunar surface and subsequently intercepted by the spacecraft as it orbits overhead. Under such 
circumstances, it may be possible to test a variety of properties and assumptions associated with various transient plume 
models or to infer certain aspects regarding lunar surface properties. 
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INTRODUCTION 

Investigators analyzing mass transport issues associated with contamination control engineering for various 
spacecraft programs sometimes encounter modeled phenomena whose behaviors are very difficult to substantiate in 
detail through ground-based testing and measurement. One such phenomenon is associated with transient 
descriptions of plume expansions from chemical rockets into high vacuum. [1-3] Very few thermal vacuum 
facilities exist that are capable of keeping high-temperature combustion or decomposition gaseous products, such as 
hydrogen, from building up during thruster testing to the point that they subsequently disturb a free expansion or 
application-specific physical features of interest. [4] 




The Lunar Atmosphere Dust Environment Explorer (LADEE) spacecraft is being designed for a mission 
featuring low altitude orbits of the Moon to take relevant ambient measurements before that environment becomes 
altered by future activities associated with human exploration (Fig. 1). Depending on orbital insertion timing, the 
observatory’s mission could last over four months, beginning in 2013. [5,6] In some cases, the orbital altitude may 
dip as low as 20 km above the lunar surface. [5] LADEE’s instrument suite includes a sensitive Neutral Mass 
Spectrometer (NMS) which will be oriented in the spacecraft’s ram direction to measure tenuous atmospheric 
constituents. Ambient density levels vary widely based on local solar illumination, but may exceed 10 5 
molecules/cm 3 during a lunar dawn. [7] NMS will be able to detect individual species concentrations at much lower 
levels, however. [8] 

LADEE’s attitude control system (ACS) includes 5-lb f bipropellant thrusters using monomethylhydrazine 
(MMH) fuel with MON-3 oxidizer, mounted in two canted pairs beneath the spacecraft’s lower deck. [5] If the ACS 
units are directed forward and toward the lunar surface, it is possible for individual species associated with their 
operations, reflected from the daytime lunar surface, to be quantitatively observed by the NMS as the observatory 
subsequently passes overhead. At lower altitudes, thruster plume fluxes to the lunar surface are maximized, as are 
reflected fluxes back to the spacecraft; also NMS may intercept peak species fluxes before disturbances from lunar 
gravity exert much influence. In this manner, NMS measurements could be used to test a variety of characteristics 
associated with transient plume models, [1,9,10] including evolution of species concentration (magnitudes, time-of- 
flight measurements, angular distribution, species separation effects), thruster exit conditions, or to infer local 
elements of lunar surface properties such as permeation or gas-surface interaction models. 


TEST CASE CONDITIONS 

LADEE will orbit the Moon in a variety of orientations and orbit altitudes. With surface values of 1737 km for 
lunar radius Ri and a gravitational acceleration constant of gi = 1.62 m/s 2 , the orbital velocity at 20 km circular 
altitude should be F 0 rb = 1 .67 km/s. [7] For this test case, a featureless strip of lunar surface was modeled along the 
ground track having width = 120 km, length = 270 km. 

For a brief thruster firing 20° below the satellite’s velocity vector with a bulk velocity of V e « 3.0 km/s, one 
would expect the average period for this pulse of gas to reach the lunar surface to be about 19.2 s in a ballistic 
trajectory due to lunar gravity, and 19.5 s in a radial trajectory without this influence. For diffusely-reflected gas of 
a dominant bipropellant thruster combustion product such as water vapor from the lunar surface at a daylight 
temperature of 380 K, [7] the time to reach 20 km based on its average thermal velocity and average direction would 
be 31.1 s with lunar gravity and 29.9 s without. 

Since ballistic perturbations from radial trajectories increase with travel time, it was decided acceptable to 
approximate conditions by neglecting the influence of gravity and limiting the period of study to one minute. Orbit 
altitude and ACS thruster firing angle were chosen to produce easily measurable peaks of dominant ACS thruster 
gases within this period. 

For the sake of simplicity, changes in LADEE’s altitude due to ACS operations were neglected. NMS will 
generally face the ram direction when collecting measurements, so only molecules reaching LADEE’s position from 
forward locations constitute the local ensemble density measured by NMS. Interactions between reflected plume 
gases and the ambient atmosphere (exosphere) were neglected since the level of rarefaction would be so 
considerable. 

LADEE’s attitude control system (ACS) includes four 5-lb f bipropellant thrusters using monomethylhydrazine 
(MMH) fuel and MON-3 oxidizer. Table 1 presents mass fractions for dominant species assumed herein. [11] 


TABLE 1. Dominant ACS thruster plume species mass fractions. [1 1] 


Species 

n 2 

h 2 o 

CO 

co 2 

h 2 

Mass Fraction 

0.426 

0.294 

0.178 

0.086 

0.016 


For assumed performance parameters and frozen flow properties, exit conditions featured a bulk velocity V e ~ 3.0 
km/s and a static temperature T e « 550 K. In the test case studied, it was assumed two forward-facing ACS units 
would fire for one second at angles 20° below LADEE’s velocity vector towards the lunar surface. 




MODEL FORMULATION 


In order to determine the magnitude of impact lunar surface-reflected ACS plume gases could have on NMS 
measurements, considering the vast scales of distances and levels of rarefaction involved, it was decided to use a set 
of point source, free molecule (FM) equations to approximately describe both the plume gas approaching the lunar 
surface as well as that diffusely reflected from it. [10] To develop this model, the collisionless Boltzmann equation 
had been solved for a directionally-constrained point source Q x meant to describe directed flow from a nozzle exit 
over 2^-steradians centered on the source normal. 
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FIGURE 2. Schematic representation of various quantities and angles used in analytical model. 

Q i represents directed flow from a Lambertian, thermal velocity distribution superimposed on convective exit 
velocity u e for mass flow rate m , (3 = 1 J4lRT , and speed ratio s = /? w e . Normal n represents the orientation of a 
local starting surface element, and v- n emphasizes the imposed directional constraint. Generally u e is not aligned 
with h , the angle between the two being defined by (j> Q (Fig. 2). Local angle 0 is measured between variable 
position x (distance r, experiencing local velocity v) and n , and angle 0 is measured between u e and x. 

Another FM source had been developed by Narasimha in conjunction with an analytical study of flow through a 
circular orifice into high vacuum. [9] This source, denoted by g N , is derived from assuming the free expansion of 
a suddenly unbounded gas: 

0N = f= S (*)»* W ex p(“ P 2 ( v - “e f )• ( 2 ) 


Clearly Q\ differs from g N by a directional constraint that produces a distinctly different character in the limit 
of s = 0; expressions developed using Q x are consistent with outgassing or surface desorption, while those using g N 
lack directional dependence. Although a recent study demonstrated the consequences of substituting g N for Q\ in 
analyzing molecular beam and pulsed laser ablation data, [10] it appeared possible to test the relative fidelity of 
these models using the NMS measurement scenario described above. 

Earlier experience indicated that a single-second thruster pulse evolution could be replaced by an instantaneous 
pulse source. [1,2] When m in Eqns. (1) & (2) is described by a Dirac Delta function Am 5 (f) , one may solve for 

successive moments of / to find pulse-mode response expressions for quantities like density p and mass flux 6 . 
For fixed geometries with a target surface element at location x relative to the plume source, incident fluxes make 
angle y/ with the local receiver normal, as shown in Fig. 2 above. 
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In Eqns. (3a) and (3b), z= a - w, a = fir /t , and w = s cos 0 . For Eq. (4), one may substitute either form for 
density to obtain the corresponding mass flux. An impermeable, inert lunar surface was assumed for this test case, 
therefore species mass flux 6, reflected from a local surface element with area A A at any given time t, was equal to 
that reaching it from the incident plume. The surface element center then became a new source using that incident 
flux in either Eq. (3 a) or (3b) with magnitude Am * 6(f) AA At to create a contribution to the reflected density at all 
later timesteps measured at NMS moving with V orb . 

If desired, it would be possible to model the effect of a permeable lunar regolith and/or the effects of surface 
chemistry by replacing the mass conservation statement with suitable mathematical models. Other gas-surface 
interaction models could be substituted for Eqns. (3a) and (3b) to produce scattered reflected fluxes. In particular, 
Eq. (3a) is consistent with Maxwell’s model for diffuse scattering. [3] 


RESULTS 

The test case was run using both Q { and Q^type models in pulse mode. Results include logarithmic contour 
maps of overall mass density across a featureless section of the lunar surface from ACS operations, as well as 
transient plots of surface-reflected number density levels intercepted by a forward-facing NMS instrument. The 
contour maps take advantage of symmetry about the axis of LADEE’s ground track. 



FIGURE 4. Test case mass density impingement contour maps at lunar surface after gn pulse: (a) 20 s, (b) 40 s, & (c) 60 s. 

In Fig. 3, a sequence of snapshots conveys the transient development of density across a featureless lunar surface 
due to a one-second LADEE ACS thruster pulse in 20-second increments. A maximum peak value of 2.7x1 O' 15 


g/cm 3 was reached along the centerline at 15 s. Figure 3 also illustrates how the plume convected along LADEE’s 
orbital direction as its footprint spread with time. 

Each thruster was modeled using a single Q x source. Since these only provide non-zero values over 2n 
steradians centered on the thruster nozzle axis, there is an abrupt cutoff in the contours that propagates at V orb due to 
this simplification. A more realistic analysis would incorporate locally varying conditions across a freezing surface 
capturing the thruster near-field expansion into the backflow regime. That data would be used to create a network of 
point sources to describe the expansion and create a better representation of the surface density map in that region. 

Solutions for species mass fluxes were also computed to propagate their subsequent presence at the LADEE 
NMS aperture. They exhibit behavior similar to the evolution of density and are not presented here. 

Figure 4 presents results using the Q N source, which share surprising similarities to Q x results depicted in Fig. 3. 
Near the maximum flux location at the time of peak intensity for the source, the maximum density level was 96 

percent as high as that for Q x results. However, it does not exhibit the cutoff for positions beyond right angles off 
the thruster plume axis. Since these are molecular flow models, one should wonder how a neutral gas flow directed 
forward via a nozzle can find itself at such high angles in the absence of any physical mechanism. 

It was noticed that the elapsed times for peak individual species mass fluxes to reach the lunar surface using 
either Q x or Q N type sources occurred more rapidly than expected (19-20 s) based on the geometry and nozzle exit 
speed, and it appears the basis for this is derived from the models themselves. Reviewing Eqns. (3) & (4), noting 
that flux M varies as t D , setting the time derivative to zero for peak flux yields 
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For w = s on the plume centerline, ^maxfiux -» r/V t as 5 — » oo . For any finite value of s however, the time to 
reach maximum flux actually decreases due to the component of kinetic energy associated with the non-zero thermal 
distribution. To observe this effect, it was decided to run a case using the Q x source where the nozzle exit static 
temperature T e was decreased from 550 K to 55 K. Results are depicted in Fig. 5 below. Differences between those 
results (Fig. 5) and Fig. 3 are much more substantial than those between Figs. 3 & 4. For the lower T e case, peak 
fluxes tend to occur between 19-20 s, and flux distributions are more tightly concentrated about those peaks. 



FIGURE 5. Surface mass density for Q x pulse, same conditions as Fig. 3, except for low T Q : (a) 20 s, (b) 40 s, & (c) 60 s. 


For the final phase of each run, lunar surface fluxes for each species at every timestep and location were used to 
create transient, thermally-accommodated, diffuse (s = 0), Q x and Q N sources to calculate their contributions to 
number density measured by the LADEE NMS as it would subsequently pass overhead. Results are presented in 
Fig. 6a ( Q x vs. Q N ) and Fig. 6b (Q u T & = 550 K vs. T e = 55 K) below. 

In Fig. 6a, Qrgenerated estimates gave rise to higher levels than Q N values, with greater discrepancies occurring 
for higher molecular weight species, along with steeper gradients. The basic deficit of Qn fluxes relative to Q x is 
due to the nature of the Q N source, which is spherically-symmetric for s = 0. As a consequence, half its material is 
directed below the lunar surface, violating the conservation assumption for mass fluxes used in this study. Some 
investigators using g N -type sources constrain the distribution function to exist only in the outbound hemisphere, 
[12] but this fix does not produce a Lambertian cosine distribution. Nevertheless, it could be difficult to distinguish 
which source model exhibits closer fidelity using NMS density measurements once such data becomes available. 


In Fig. 6b, the effect of higher speed ratio and fixed V e is readily apparent. The delay in peak surface fluxes 
translates to a longer period before reaching maximum density levels, and the tighter footprint for lunar surface mass 
flux produces higher peak density levels once they do occur. 



TIME ELAPSED FROM THRUSTER FIRING [s] 

(a) 



TIME ELAPSED FROM THRUSTER FIRING [s] 

(b) 


FIGURE 6. Comparison of predicted species contributions measured by LADEE NMS. (a) Q x (solid) vs. Q N (dashed), 
(b) effect of 7 e on Q x ; baseline ( T c = 550 K, solid) vs. 10x lower value ( T c = 55 K, dashed). 


CONCLUDING REMARKS 

A study was performed to test the possibility that LADEE ACS thruster operations could be observed by its 
NMS instrument after reflection from the lunar surface. Results indicate such measurements are possible, and the 
study explored the possibility that the fidelity of various transient source flow plume models could be rated. It may 
also be possible to similar assessments for models describing the effects of a permeable lunar regolith, surface 
chemistry, and gas-surface interaction models. 
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